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Abstract 

A  numerical  model  of  a  three-stage  magnetic  modulator  for  cop¬ 
per  lasers  is  described.  The  model  follows  in  detail  the  hysteresis 
behavior  of  the  magnetic  switches  and  a  transformer.  The  laser  is 
treated  via  a  simple  model  of  field  diffusion.  The  calculated  voltage 
across  the  laser  and  current  through  the  laser  are  shown  to  compare 
favorably  with  experiment. 


Introduction 


This  paper  presents  simulations  of  a  magnetic  modulator  which 
was  developed  to  drive  copper  lasers  used  in  isotope  separation  [1]. 
The  circuit  is  shown  schematically  in  Fig.  1;  parameter  values  are 
listed  in  Table  1.  The  modulator  uses  a  thyratron  or  SCR  stack  and 
three  stages  of  magnetic  compression  to  provide  a  40-45  kV  pulse 
to  the  laser  at  a  repetition  rate  of  4.5  kHz.  It  operates  at  an  input 
power  of  28-36  kW.  Each  stage  has  a  saturable  inductor  consisting 
of  several  cores  made  of  the  ferrite  C/7D.  There  is  also  a  1:4  step-up 
transformer  wound  with  Metglas  2605S3A. 

Initially  the  capacitor  Co  is  charged  to  a  typical  voltage  of  18  kV. 
The  magnetic  switches  are  held  in  negative  satutation  by  reset  cur¬ 
rents.  All  other  voltages  and  currents  are  zero.  After  the  thyratron 
or  SCR  stack  closes,  Co  discharges  upon  Ci  with  a  charge  time  of 
approximately  900  ns.  During  this  time  the  first  magnetic  switch 
holds  off  the  current  by  presenting  a  large  inductance  to  the  circuit. 
It  then  saturates,  allowing  the  energy  to  flow  to  C2,  which  has  a 
charge  time  of  350  11s.  When  the  second  switch  saturates,  the  energy 
flows  to  C3  (charge  time  of  90  ns),  and  the  voltage  is  stepped  up 
by  a  factor  of  4  by  the  transformer.  After  saturation  of  the  third 
switch,  the  voltage  is  transferred  to  the  peaking  capacitor  (rise  time 
of  30  ns)  and  then  to  the  laser  load.  It  then  reflects  back  through 
the  circuit,  undergoing  expansion  in  successive  stages. 

The  reset  circuit,  not  shown  in  Fig.  1,  consists  of  a  rod  passing 
through  each  core,  connected  in  series  with  a  fixed  voltage  source,  a 
large  inductor,  and  a  resistor. 

The  task  of  the  model  is  to  describe  the  detailed  behavior  of  the 
voltages,  currents,  and  fields  in  the  modulator  system.  We  shall  see 
that  this  can  be  done  in  quite  reasonable  accordance  with  experi¬ 
ment. 


Magnetic  Model 

We  employ  a  simple  core  model  [2]  in  which  the  magnetic  field 
within  the  cores,  in  both  the  switches  and  the  transformer,  is  taken 
to  be  azimuthal  and  to  depend  only  on  time,  and  the  displacement 
current  is  ignored.  The  voltage  across  each  core  then  satisfies  an 
equation  of  the  form 


V(t)  =  L(t) 


dl(t) 

dt 


(1) 


where  the  time-dependent  inductance  is  directly  proportional  to  the 
permeability  /r(i)  within  the  core.  We  neglect  the  radial  dependence 
of  the  fields,  which  is  an  effect  of  order  10%.  The  magnetic  field  itself 
evolves  according  to 


m  =  ^(m),  (2) 

with  Nw  the  number  of  windings  and  f  an  average  radius.  Further 
details  of  this  model  are  given  in  the  Appendix. 


Fig.  1.  Schematic  modulator  circuit.  The  reset  circuit  is  not  shown 

Table  1.  Modulator  parameter  values.  Resistances  listed  next  to 
capacitances  are  nominal  ESR  values. 


Stage  1 

Stage  2 

Stage  3 

Transformer 

Material 

C/7D 

C/7D 

C/7D 

2605S3A 

Number  of  cores 

10 

8 

10 

1 

Number  of  windings 

5 

2 

2 

4  (secondary) 

Inner  radius  (cm) 

4.57 

4.57 

4.57 

4.45 

Outer  radius  (cm) 

7.37 

7.37 

7.37 

7.65 

Core  height  (cm) 

1.27 

1.27 

1.27 

10.16 

Saturated  L  (nH) 

637 

67.8 

450 

11  (primary) 

Packing  factor 

.95 

.95 

.95 

.65 

C0  =  40.5  nF  RCo 

Ci  =  40.5  nF  Rci 

Ci  =  40.5  11F  Rc2 

=  67  mfi 
=  67  mfi 
=  67  mfi 

Lo  = 

5  /iH 

C3  =  2.5  nF  Rc3  =  400  mfi 

Reset  circuit  parameters:  Vre  =  12  v,  Rre  =  1  fl,  Lre  =  1.2  mH 
Laser  parameters:  a  —  4  cm,  l  =  273  cm,  Lc  =  400  nH 
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Hysteresis  properties  are  contained  in  the  permeability  /r(t).  In 
the  model,  we  employ  a  hysteresis  theory  developed  by  Hodgdon  [3], 
according  to  which  the  permeability  varies  with  B,  B ,  and  H  in  the 
maimer  summarized  in  the  Appendix.  This  theory  reproduces  many 
of  the  known  features  of  hysteresis,  such  as  loop  closure  near  satu¬ 
ration,  the  formation  of  minor  loops,  and  loop  broadening  with  in¬ 
creasing  B.  The  theory  requires  several  material  parameters;  those 
appropriate  to  C/7D  and  Metglas  2605S3A  are  given  in  the  Ap¬ 
pendix. 

In  previous  applications  [2,4],  this  combination  of  core  model  and 
hysteresis  model  was  shown  to  agree  satisfactorily  with  experiment. 

We  employ  a  simple  model  of  the  trigger  switch  (thyratron  or  SCR 
stack).  For  positive  currents,  it  is  treated  as  a  negative  500  v  source, 
to  account  for  the  forward  drop  across  the  switch.  For  negative 
currents,  it  is  modeled  as  a  400  ohm  resistor  in  series  with  the  same 
source. 


where  the  value  of  the  electric  field  in  the  outermost  cell  is 

£/v-i  =  7(4-iSn  —  i).  (8) 

2  Poo 

Here  Vpk  is  the  voltage  across  Cpk  and  its  equivalent  series  resistance, 
and  l  is  the  length  of  the  laser. 

The  final  model,  therefore,  consists  of  a  system  of  15  circuit  vari¬ 
ables  (5  voltages,  6  currents,  and  4  core  magnetic  fields),  along  with 
the  JV  —  1  =  8  magnetic  field  variables  in  the  laser.  Because  of  the 
core  magnetic  model,  this  is  a  stiff  system  of  ODEs.  We  advance  it 
in  time  via  the  solver  LSODKR,  a  variant  of  the  package  LSODE  [5], 
which  employs  Krylov  methods  in  numerical  iterations.  The  solver 
also  has  root-finding  capabilities  which  are  used  in  transitions  be¬ 
tween  branches  of  the  hysteresis  curves. 

Results 


Laser  Model 


From  the  point  of  view  of  the  modulator,  the  laser  plays  the  role 
of  a  complicated  load.  For  purposes  of  understanding  modulator 
behavior,  it  suffices  to  treat  the  laser  as  a  cylindrical  plasma  with 
constant  resistivity  p,  in  series  with  the  fixed  inductance  Lc  of  the 
laser  container.  Obviously  this  model  of  the  laser  cannot  predict 
detailed  properties  such  as  the  optical  power. 

We  adopt  a  one-dimensional  model  in  which  the  electric  field 
E(r,t)  points  in  the  z  direction  and  the  magnetic  field  B(r,t)  is 
purely  azimuthal.  In  the  absence  of  displacement  current,  the  mag¬ 
netic  field  then  satisfies  a  diffusion  equation  of  the  form 


dB_  _  d  p  d(rB) 
dt  dr  fj.0r  dr  ’ 


with  a  boundary  value  at  the  radius  a  corresponding  to  the  plasma 
current:  B(a,t)  =  fi0Ip(t)/2ira.  At  each  instant,  the  electric  field  is 
given  as 


P  _  P  d(rB) 
H  or  or 


(4) 


and  thus  is  not  an  independent  dynamical  quantity.  The  appropriate 
value  of  the  resistivity,  obtained  by  comparing  a  more  detailed  model 
with  experiment,  turns  out  to  be  approximately  .3  ohm  cm. 


Method  of  Solution 


It  is  convenient  to  formulate  the  model  as  a  system  of  ordinary 
differential  equations  suitable  for  an  ODE  solver.  The  circuit  equa¬ 
tions  are  already  in  this  form,  although  the  hysteresis  model  presents 
extreme  nonlinearities  and  jumps  in  the  form  of  the  equations. 

To  handle  field  diffusion  in  the  plasma,  we  follow  a  conventional 
approach  in  which  the  magnetic  field  is  assigned  values  at  N  nodal 
points  distributed  evenly  from  the  center  to  the  edge.  In  practice,  it 
was  found  sufficient  to  take  N  =  9.  At  points  2  through  N  —  1,  the 
field  then  satisfies  the  three-point  difference  equation 

Bi  =  [4d"  5;+ 1  +  a“  Bi  +  ai  f?;_ ij ,  (5) 


with 


2(»-l)2 

(<-*)(<-!)’ 


i-  2 

M' 


(6) 


and  (S  the  spacing  between  points.  Of  course  B i  is  identically  zero. 
A  consistent  approximation  for  the  edge  value,  obtained  in  part  by 
differentiating  the  boundary  condition,  is 


Bn  =  {Vpk  -  £EN_i)/(2iraLc/ fj.Q  +  S 1/ 2),  (7) 


After  the  voltage  pulse  is  compressed  and  reflects  back  to  Co,  a 
recovery  system  stores  energy  for  the  next  pulse.  Since  this  system  is 
not  treated  by  the  model,  simulations  are  stopped  when  the  voltage 
across  Co  reaches  a  minimum.  This  occurs  about  2.6  /is  after  the 
thyratron/SCR  is  fired. 

Figure  2  shows  calculated  voltage  waveforms  for  Co,  Ci,  and  C2 
during  this  interval.  As  one  would  expect,  the  900  ns  falloff  in  voltage 
across  Co  is  matched  very  closely  by  the  rise  across  Ci .  After  the 
first  magnetic  switch  fires,  the  fall  across  Ci  (350  ns)  is  matched 
by  the  rise  across  C'2.  After  the  second  switch  fires,  the  voltage  is 
transferred  from  C2  to  C3  in  90  11s.  Note  that,  after  compression,  the 
waveforms  expand  as  the  energy  reflects  back  to  Co-  These  features 
are  in  excellent  agreement  with  expectations  based  on  the  design  of 
the  compression  system. 

The  calculated  hysteresis  behavior  of  the  three  magnetic  switches 
is  shown  in  Figures  3,  4,  and  5.  Since  successive  stages  experience 
increasingly  higher  field  frequencies,  the  B-H  loops  broaden  in  width 
and  develop  minor  loops.  The  widths,  as  measured  at  two-thirds  the 
saturation  value  of  B,  are  160  At/m,  260  At/m,  and  520  At/m,in 
the  respective  stages.  All  three  stages  saturate  twice  in  the  positive 
sense.  In  stage  1  the  second  saturation  occurs  outside  the  plotted 
domain,  but  in  the  other  two  stages  the  loop  corresponding  to  the 
second  saturation  is  visible.  The  dissipated  hysteresis  powers  are 
175  W,  270  W,  and  710  W,  respectively. 

The  voltage  across  the  laser  is  depicted  in  Fig.  6.  Note  that  it 
peaks  at  37  kV,  drops  within  90  ns  to  a  minimum  of  -37  kV,  and 
then  slowly  rings  down  with  a  period  of  about  240  ns.  This  behavior 
agrees  quite  well  with  the  measured  voltage  trace  [6],  except  that 


305 


the  latter  shows  somewhat  less  dissipation.  The  measurement  also 
reveals  an  initial  leakage,  which  could  be  modeled,  if  desired,  by 
placing  a  capacitor  across  the  third  stage  switch. 

Figure  7  shows  the  current  through  the  laser.  The  calculated  cur¬ 
rent  peaks  at  37  kA  and  then  rings  down.  The  measured  waveform 
reaches  a  maximum  of  about  60%  of  the  calculated  value,  but  there 
may  be  a  substantial  experimental  error  in  this  region.  Model  and 
experiment  are  reasonably  close  in  the  subsequent  decay. 

Finally,  Figure  8  shows  the  calculated  electric  field  in  the  laser  as 
a  function  of  time,  at  three  radial  locations.  The  field  is  largest  near 
the  edge,  reaching  a  maximum  of  38  v/cm,  then  decays  as  it  diffuses 
inward. 


Conclusions 


These  simulations  show  that  the  modulator  system  can  be  modeled 
with  quite  reasonable  accuracy.  The  code  has  served  as  a  useful 
guide  in  understanding  circuit  behavior  and  in  the  design  of  possible 
modifications  to  the  device. 


Fig.  5.  Hysteresis  loop  of  third  magnetic  switch. 


Fig.  3.  Hysteresis  loop  of  first  magnetic  switch. 


Fig.  6.  Voltage  across  laser.  The  dotted  line  indicates  Chang’s  data. 


Time  (microsec) 


Fig.  4.  Hysteresis  loop  of  second  magnetic  switch. 


Fig.  7.  Laser  current.  The  dotted  line  indicates  Chang’s  data. 


306 


of  B  leads  to  hysteresis.  The  material  functions  /  and  g  have  the 
form 


Time  (microsec) 


Fig.  8.  Electric,  field  in  laser,  at  three  radii. 
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Appendix:  Details  of  Core  Magnetic  Model 

Neglecting  the  radial  dependence  of  the  fields,  the  time-dependent 
inductance  of  a  core  is  [2] 


m  «  (Al) 

2ir  Tin 

where  h  is  the  height  of  a  core,  Nw  is  the  number  of  windings,  and 
and  rout  are  the  inner  radius  and  outer  radius,  respectively.  The  time 
derivative  of  the  magnetic  field,  Eq.  (2),  is  obtained  by  combining  the 
constitutive  relation  6B  =  f.i6H  with  Ampere’s  law.  The  appropriate 
average  radius  is  f  =  (rou(  -  rin)/\n(rout/rin). 

For  convenience,  we  collect  here  the  needed  elements  of  the  Hodg- 
don  theory  of  hysteresis  [3],  The  permeability  varies  with  B ,  B,  and 
H  according  to 

KB,  H,  B)  =  {a  sign(B)  [/(B)  -  H]  +  g(B,  B)f\  (A2) 

where  a  is  a  constant  and  /  and  g  are  functions  particular  to  the 
magnetic  material.  The  dependence  of  the  permeability  on  the  sign 


m  = 


9(B,B) 


The  rate 


(  A1[tan(AaB)r  +  UB, 

<  Ai[tan(A2B(,)]J’  +  fexBb  +  {B  -  Bb)/fis, 

{ -f(-B), 


f'(B), 


0  <B<Bb] 
B  >  Bb ; 

B  <  0; 

(A3) 

if  |B|  <  Bci; 

otherwise. 

(A4) 


dependence  is  governed  by  c(B),  which  is  given  by 


fl+<h|£|,  |B|<Bi; 

C(B)  =  <  1  +  ciBi  -f  cj(|B|  —  Bi),  B,<|B|<B2; 

[  1  +  ci Bi  +  c2(B2  -  Bi)  +  c3(|B|  -  B2),  otherwise. 

(AS) 

The  values  of  the  constants  for  C/7D  are  (inks  units): 


Bc,  =  .25 
Bb  =  .2596 

fex  —  9 

P  =  1 


Aj  =  39.36 
A2  =  5.984 
A3  =  -.5579 
A4  =  .01043 


cj  =  1.1  x  10-6 
c2  =  4.8  x  10~6 
c3  =  ... 

~  fi0 


Bj  =  4.2  x  10s 
B2  =  oo 
B3  =  ... 
a  =  10 


in  which  c3  and  B3  are  not  needed  since  only  the  first  two  branches 
of  c(B)  are  used.  For  Metglas  2605S3A,  the  constants  are: 


Bci  =  1-4 

A,  =  .08734 

ci  =  4.455  X  10~4 

Bj  =  1.1  x  105 

Bb  =  1.45 

A2  =  1.0778 

c2  =  7.868  X  lO"5 

B2  =  3.3  X  106 

Sex  =  -7958 

A3  =  -8.944 

c3  =  5.833  x  10"6 

B3  =  8.7  x  106 

P  =  2 

A4  =  .2770 

/is  =  .9324/iq 

a  =  10 

References 

[1]  E.  G.  Cook  et  al.,  “High  Average  Power  Magnetic  Modulator  for 
Copper  Lasers,”  Proc.  Eighth  IEEE  International  Pulsed  Power 
Conference,  San  Diego,  CA,  June  1991,  p.  537. 

[2]  C.  D.  Boley  and  M.  L.  Hodgdon,  “Model  and  Simulations  of 
Hysteresis  in  Magnetic  Cores,”  IEEE  Trans.  Magnetics  25,  3922 
(1989). 

[3]  M.  L.  Hodgdon,  “Applications  of  a  Theory  of  Ferromagnetic  Hys¬ 
teresis,”  IEEE  Trans.  Magnetics  24,  218  (1988);  M.  L.  Hodgdon, 
“Mathematical  Theory  and  Calculations  of  Magnetic  Hysteresis 
Curves,”  IEEE  Trans.  Magnetics  24,  3120  (1988). 

[4]  C.  D.  Boley,  “Simulations  of  a  Magnetic  Pulse  Compressor,” 
Proc.  Eighth  IEEE  International  Pulsed  Power  Conf.,  San  Diego, 
CA,  June  1991,  p.  704. 

[5]  A.  C.  Ilindmarsh,  “ODEPACK,  a  Systematized  Collection  of 
ODE  Solvers,”  in  Scientific.  Computing,  R.  S.  Stepleman  et  al. 
(eds.),  North-Holland,  Amsterdatn,  1983  (Vol.  I  of  IMACS  Trans¬ 
actions  on  Scientific  Computing),  pp.  55-64. 

[6]  J.  J.  Chang,  private  communication. 


307 


